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A Possible Implementation of a Direct Coupling Coherent Quantum 

Observer 


Ian R. Petersen and Elanor H. Huntington 


Abstract — This paper considers the problem of implementing 
a previously proposed direct conpling qnantnm observer for a 
closed linear quantum system. This observer is shown to be able 
to estimate some bnt not all of the plant variables in a time 
averaged sense. The paper proposes a possible experimental 
implementation of the observer plant system using a non¬ 
degenerate parametric amplifier. 

I. Introduction 

A number of papers have recently considered the problem 
of constructing a coherent quantum observer for a quan¬ 
tum system; see [l]-[3]. In the coherent quantum observer 
problem, a quantum plant is coupled to a quantum observer 
which is also a quantum system. The quantum observer is 
constructed to be a physically realizable quantum system 
so that the system variables of the quantum observer con¬ 
verge in some suitable sense to the system variables of the 
quantum plant. The papers [4]-[7] considered the problem of 
constructing a direct coupling quantum observer for a given 
quantum system. 

In the papers [1], [2], [4], [6], the quantum plant under 
consideration is a linear quantum system. In recent years, 
there has been considerable interest in the modeling and 
feedback control of linear quantum systems; e.g., see [8]- 
[10]. Such linear quantum systems commonly arise in the 
area of quantum optics; e.g., see [11], [12]. For such linear 
quantum system models, an important class of quantum 
control problems are referred to as coherent quantum feed¬ 
back control problems; e.g., see [8], [9], [13]-[18]. In these 
coherent quantum feedback control problems, both the plant 
and the controller are quantum systems and the controller is 
typically to be designed to optimize some performance index. 
The coherent quantum observer problem can be regarded 
as a special case of the coherent quantum feedback control 
problem in which the objective of the observer is to estimate 
the system variables of the quantum plant. 

In this paper, we consider the situation as in papers [4]-[7] 
in which there is only direct coupling between quantum plant 
and the quantum observer. In these papers, both the quantum 
plant and the quantum observer are assumed to be closed 
quantum systems which means that they are not subject to 
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quantum noise and are purely deterministic systems. This 
leads to an observer structure of the form shown in Figure 
[T] In these papers, it is shown that a quantum observer can 
be constructed to estimate some but not all of the system 
variables of the quantum plant. Also, the observer variables 
converge to the plant variables in a time averaged sense rather 
than a quantum expectation sense such as considered in the 
papers [1], [2]. 



Fig. 1. Coherent Observer Structure with Direct Coupling. 

In this paper, we concentrate on the result presented in 
[4] for the case in which the quantum plant is a single 
quantum harmonic oscillator and the quantum observer is 
a single quantum harmonic oscillator. For this case, we 
show that a possible experimental implementation of the 
augmented quantum plant and quantum observer system may 
be constructed using a non-degenerate parametric amplifier 
(NDPA) which is coupled to a beamsplitter by suitable choice 
of the amplifier and beamsplitter parameters. 

II. Quantum Linear Systems 

In this section, we describe the class of closed linear 
quantum systems under consideration; see also [4], [6], [8], 
[15], [19]. We consider linear non-commutative systems of 
the form 

x{t) = Ax{t)-, a;(0) = xq (1) 

where A is a real matrix in and x{t) = 

[ Xi{t) ... Xn{t) Y' is ^ vector of self-adjoint possibly 
non-commutative system variables; e.g., see [8]. Here n is 
assumed to be an even number and ^ is the number of modes 
in the quantum system. 

The initial system variables a;(0) = xq are assumed to 
satisfy the commutation relations 

[a;j(0),a;fe(0)] = 2i0jfc, j, fc = l,...,n, (2) 

where 0 is a real antisymmetric matrix with components 
Qjk. Here, the commutator is defined by [A, B] = AB—BA. 
In the case of a single degree of freedom quantum particle, 
x = {xi,X 2 )'^ where xi = q is the position operator, 
and X 2 = p is the momentum operator. The commutation 
relations are \q,p] = 2i. Here, the matrix 0 is assumed to 





be of the form 0 = diag( J, J,... ,J) where J denotes the 
real skew-symmetric 2x2 matrix 


A linear quantum system O is said to be physically 
realizable if it ensures the preservation of the canonical 
commutation relations (CCRs); 

x(t)x(tY' — {x{t)x{t)'^Y' = 2iQ for all f > 0. 

This holds when the system O corresponds to a collection of 
closed quantum harmonic oscillators', see [8]. Such quantum 
harmonic oscillators are described by a quadratic Hamilto¬ 
nian H = ^x"^Rx, where i? is a real symmetric matrix. 

In the proposed direct coupling coherent quantum ob¬ 
server, the quantum plant is a single quantum harmonic 
oscillator which is a linear quantum system of the form ([T]i 
described by the non-commutative differential equation 


ipil') — — ^Op: 

Zp{t) = CpXp{t) (3) 


where Zp{t) denotes the vector of system variables to be 
estimated by the observer and Ap G 


2x2 


Cp G 


plx2 


It 


is assumed that this quantum plant corresponds to a plant 


Hamiltonian Up = ^x'ERpXp. Here Xp = 


Qp 

Pp 


where Qp 


is the plant position operator and pp is the plant momentum 
operator. 

We now describe the linear quantum system of the form 
O which will correspond to the quantum observer; see 
also [8], [15], [19]. This system is described by a non- 
commutative differential equation of the form 


Xo{t) = AoXo{t)', Xo{0) = XQo', 

Zo{t) = CoXo{t) (4) 


HI. Direct Coupling Distributed Coherent 
Quantum Observer 

Following [4], [6], we assume that Ap = 0 in Q. This 
corresponds to Rp = 0 in the plant Hamiltonian. It follows 
from Q that the plant system variables Xp{t) will remain 
fixed if the plant is not coupled to the observer. However, 
when the plant is coupled to the quantum observer this will 
no longer be the case. In addition, we construct the observer 
as in [4] so that 

> 0 and i?c = Cj/3 (6) 

where 

13 G and CoR~^/3'^ = -1. (7) 

With this construction of the quantum observer, the fol¬ 
lowing result was established in [4]. 

Theorem 1: Consider a quantum plant of the form Q 
where Ap = 0. Then the matrices Ro > 0, Rc, Co satisfying 
the conditions (|6]l, (|7]i will define direct coupled quantum 
observer such that for the resulting augmented plant-observer 
system, the quantity Zp{t) is constant and 

lim i / {zp- Zo{t))dt = 0. (8) 

T-^oo 1 Jq 

IV. A Possible Implementation oe the Plant 
Observer System 

In this section, we describe one possible experimental 
implementation of the plant-observer system given in the 
previous section. The plant-observer system is a linear 
quantum system of the form O with Hamiltonian R — 
\x'^Rx = ^XoRqXo + XpRcXo where the conditions (|6l), 
© are satisfied. In particular, we assume Ro = 2ujoI > 0 
and hence 

n = ujox'l Xo AXpCp^Xo- ( 9 ) 

Also, the condition © becomes 


where the observer output Zo{t) is the observer estimate 


and Ao G 


p2x2 


, Co G 


plx2 


. Here Xo = 


do 

Po 


where 


Qo is the observer position operator and po is the observer 
momentum operator. We assume that the plant variables 
commute with the observer variables. The system dynamics 
@ are determined by the observer system Hamiltonian 
which is a self-adjoint operator on the underlying Hilbert 
space for the observer. For the quantum observer under 
consideration, this Hamiltonian is given by a quadratic form: 
Ro = \x'^RoXo, where Ro is areal symmetric matrix. Then, 
the corresponding matrix Ao in © is given by 


Ao = 2JRo. 


(5) 


In addition, we define a coupling Hamiltonian which de¬ 
termines the coupling between the quantum plant and the 
quantum observer: 


Cof + 2wo = 0. (10) 

In order to construct a linear quantum system with a 
Hamiltonian of this form, we consider an NDPA coupled 
to a beamsplitter as shown schematically in Figure |2j e.g., 
see [12]. 




Fig. 2. NDPA coupled to a beamsplitter. 














A linearized approximation to the NDPA is defined by a 
quadratic Hamiltonian of the form 

Ui = ^{ea*b* -e*ab)+uJob*b 

where a is the annihilation operator corresponding to the 
first mode of the NDPA and b is the annihilation operator 
corresponding to the second mode of the NDPA. These 
modes will be assumed to be of the same frequency but with 
a different polarization with a corresponding to the quantum 
plant and b corresponding to the quantum observer. Also, e 
is a complex parameter defining the level of squeezing in the 
NDPA and Wo is the detuning frequency of the b mode in the 
NDPA. The a mode in the NDPA is assumed to be tuned. 
In addition, the NDPA is defined by the vector of coupling 

. Here 7 > 0 is a scalar parameter 

y/lb \ 

determined by the reflectance of the mirrors in the NDPA. 

From the above Hamiltonian and coupling operators, we 
can calculate the following quantum stochastic differential 
equations (QSDEs) describing the NDPA: 


operators L = 


da 

db 


0 f 

f 0 

1 
2 


a 

b* 


dJ^out 

dBout 


w 2 r- 

Vi 0 

0 y 

Vt 0 
0 Vt 



dt 


• 

a 

o 


b 


dA ' 


dB 


dt 


dt - 


dA 

dB 


Cll) 


e.g., see [ 10 ]. 

We now consider the equations defining the beamsplitter 


A 

B 


COS0 e 

—e®‘^sin 0 cost 



J^out 


^out 


where d and (j) are angle parameters defining the beamsplitter; 
e.g., see [20]. This implies 


j^OUt 

QOUt 


cos 9 —e ” 
sin 9 cos 9 


sin 9 



■ A ' 


B 


Substituting this into the second equation in (ITSli . we obtain 


cos 9 —e ®‘^sin0 


■ dA ' 

sin 0 cos 0 


dB 


= v^ 


dt ' 


dA 

dB 


and hence 


cos 9—1 

—e *‘^sin 0 


■ dA ' 

= V^ 

a 

sin 0 

cos 0—1 


dB 

b 


dt. 


Substituting this into the first equation in ([T 2 ]) . we obtain 


da 

db 


2 

2 

2 

0 


2 

0 

0 

7 


dt 

a 

b 


dt 


7 


2 + lUJo 

COS 0 — 1 e^^^sin^ 
—e*‘^sin 0 cos 0—1 


2(1 — COS0) 

These QSDEs can be written in the form 



a 


b 


dt. 


da 


a 

db 

= F 

b 

da* 

a* 

_ db* _ 


b* 


dt 


where the matrix F is given by 

sin 0 


F = 


0 

'ye'^^ sin 6 

2(1—COS 9 ) 

0 


2(1—cos 9) 

-lUJo 

€* 

2 

0 


2 

0 

7 e~^^ sin 9 

2(1—cos 9) 


76 ^*^ sin 9 

2(1—cos 9) 


lUJn 


It now follows from the proof of Theorem 1 in [10] that we 
can construct a Hamiltonian for this system of the form 


B. = -[a* b* a b]M 


where the matrix M is given by 


a 

b 

a* 

b* 


M = ^ (JF-FU) 


where J = 


I 0 
0 -I 


. Then, we calculate 


M = 


0 

76 *'^ sin 9 

1 —cos 9 

0 

-e* 


76 


’ sin 9 


1—cos 9 
— 2iU!o 
-e* 

0 


e 

0 

■ye"' 




1 —cos 9 


1—cos 9 
—2tuJo 


We now wish to calculate the Hamiltonian H in terms of 
the quadrature variables defined such that 


a 


Qp 

b 

= $ 

Pp 

a* 


do 

. b* _ 


. Po . 


We now assume that cos0 jtz l. It follows that we can write where the matrix $ is given by 


dA 

dB 




2(1 — COS0) 


cos 0 — 1 
sin0 


e-*'^sin 0 
cos 0 — 1 



a 


b 


$ = 


dt. 


1 i 0 0 

0 0 It 

1 -i 0 0 

0 0 1 


—I 

































































































Then we calculate 


40 


= 


T T 

Xp Xg 


R 


where the matrix R is given by 
R = 


Rc = 


and a = 

1 —COS & 


0 Rc 

R^ 2uj„I 


-3(e)-3(a) 5i(e) + K(a) 

3(e) — 3(a) 


5R(e) - 5ft(a) 
( Hence, 






Qp 

30 - 

Po ] R 

Pp 

20 - 


Qo 



. Po . 

10 - \ 

Xp 



Xo 


sr 0 


Fig. 3. Plot of the function f{9). 


R = LOoX^Xo + XpRcXo- 

Comparing this with equation (|9l), we require that 

—3(e) — 3(a) 5R(e) + 5R(a) 

5R(e) — 5ft(a) 3(e) — 3(a) 

and the conditions ©, © to be satisfied in order for the 
system shown in Figure [2] to provide an implementation of 
the augmented plant-observer system. 

We first observe that the matrix on the right hand side of 
equation (fTSli is a rank one matrix and hence, we require that 

= |a|2-|ep = 0. 

That is, we require that 


—3(e) — 3(a) 5R(e) + 5ft(a) 

5i(e) — 5i(a) 3(e) — 3(a) 



sin 0 
1 — cos 6 


e ■ 


Note that the function , ^ „ takes on all values in (—oo, oo) 

for 9 € (0, 27r) and hence, this condition can always be 
satisfied for a suitable choice of 9. This can be seen in Figure 
[3 which shows a plot of the function f{9) = 

Furthermore, we will assume without loss of generality 
that 9 G (0, tt) and hence we obtain our first design equation 


sin 0 |e| 

1 — cos 9 7 


(13) 


In practice, this ratio would be chosen in the range of G 
( 0 , 0 . 6 ) in order to ensure that the linearized model which is 
being used is valid. 

Our second design equation is obtained from (fT3 l. First we 
write Cp = [ Cpi Cp 2 ] and define the complex number 
c = Cpi+iCp 2 . The argument of this complex number arg(c) 
determines the quadrature of interest in the plant. It is straight 
forward to verify that the condition (fT2l l will be satisfied for 
some non-zero vector /3 if and only if 


arg (z(e — a*)) = arg(c). 


That is. 


arg - 

This is equivalent to 



76 sin 9 
1 — cos 9 

g-i<^ sin 0 
1 — cos 9 


^ = arg(c) 
^ = arg(c) 


TT 

2 ' 


TT 

2 ' 


(14) 


Then we write e = |e|e®^. It follows from (fTsT l that (fl4l i can 
be re-written as 


arg (e*’^ - e *‘^) = arg(c) - (15) 

This is our second design equation. 

If the design equations (fTST l and (fTSl) are satisfied then 
there will exist a non-zero vector /3 such that (fT2]) will 
be satisfied. Then, we can always find a non-zero vector 
Co such that (fTOl i is satisfied. Thus, using Theorem [T] we 
can conclude that if the the design equations ( fOl l. ( fTSl ) are 
satisfied then the corresponding direct coupled observer will 
have the desired properties. However, since the proposed 
experimental implementation of the plant-observer system 
is a closed quantum system, there is no measurement which 
can be made on the system to verify the performance of this 
system. Future research will be directed towards extending 
the theory developed in [4] to allow for open quantum 
systems in which a measurement can be made to verify the 
behaviour of the direct coupled quantum observer. 

Example 

We now illustrate the above design principles with an 
example corresponding to typical laboratory values. In this 
example, we let Cp =[ 1 0 ] corresponds to the case in 

which the position quadrature of the plant is of interest. Also, 
we choose 7 = 10® rad/s and ojo = 10® rad/s. In addition, 
we choose 



7 


which according to O corresponds to a value of 9 = 
168.6°. 

























In this case c = 1 which is purely real and the condition 
(ITSl l reduces to 



That is 

COs('!/') — COs((/)) = 0 

and 

sin('0) + sin((^) < 0. 

To satisfy these conditions, we choose ij} = (j) = —90°. Then, 
we obtain 

e = = -zlO^ 

and 

a = 10^e"*-5 = -ilO^ 

Then, it follows from (fT2T i that 

■ 2 X 10^ 0 

0 0 

Hence 

/3 = [ 2 X 10^ 0 ] . 

Therefore, if we choose Co = [ — 10 0 ], it follows 
that the condition Coll will be satisfied. Thus, with these 
parameter values, the proposed implementation will satisfy 
the conditions of Theorem [1] for a direct coupled quantum 
observer. 

V. Conclusion 

We have shown that the direct coupling quantum observer 
proposed in [4] could be at least in theory implemented 
experimentally. However, such an experiment could not 
provide experimental verification that the properties of such 
a quantum observer described in Theorem [T] are satisfied. 
In order to address this issue, future research will extend 
the results of [4] to allow for a small probe field. Then the 
theory developed in this paper will be extended to allow for 
this case. 

Another area of possible future research would be to 
analyse the performance of the proposed implementation 
of the plant-observer system without making a linearization 
assumption in the model of the NDPA. 
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